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Abstract 

The modelling of DNA mechanics under external constraints is dis- 
cussed. Two analytical models are widely known, but disagree for instance 
on the value of the torsional modulus. The origin of this embarassing sit- 
uation is located in the concept of writhe. This letter presents a unified 
model for DNA establishing a relation between the different approaches. 
I show that the writhe created by the loops of DNA is at the origin of the 
discrepancy. To take this into account, I propose a new treatment of loop 
statistics based on numerical simulations using the most general formula 
for the writhe, and on analytic calculations with only one fit parameter. 
One can then compute the value of the torsional modulus of DNA without 
the need of any cut-off. 



1 General motivations 

New experimental techniques in single molecule manipulation of DNA and su- 
percoiling control have stimulated improvements in the understanding of DNA 
mechanics ^ |21 El El ■ Surprisingly, the measurements have to be interpreted 
through a rather sophisticated model in order to extract physical constants [510 
[7]. Different approaches lead to disagreeing values of, in particular, the torsional 
modulus C of the molecule along its axis. The method used by Moroz and Nel- 
son [HI leads to the value C/ksT — 109 nm, while using the same experimental 
data the model provided by Mezard and Bouchiat [7| gives 84± 10 nm. Recently, 
another experiment performed two direct measurements with a weighted av- 
erage of 102 ± 6 nm. In this letter, I will establish the domains of validity of the 
two mentioned theoretical approaches, clarify the origin of their disagreement 
and compute a value of the torsional constant C with a new model. 

This letter is organized as follows: It starts with a short introduction to two 
widely used models for elasticity of a polymer. In order to take into account 
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DNA resistance to torsion, one of them, worm-like chain model, is improved 
with a new elastic energy term. As these models only contain a local descrip- 
tion of the molecule, the following paragraph introduces the global geometrical 
description of a single DNA molecule. Its three dimensional coiling is described 
by a quantity called writhe, whose fluctuations reduce the effective torsional 
modulus. I relate the models of Ref. [HI and [7] and study their domains of 
validity. The writhe fluctuations are interpreted in the following paragraph by 
considering the loops formed by the molecule. The last paragraph describes both 
a full numerical computation of the writhe fluctuations and a one parameter fit 
for the interpretation in terms of loops. The different theoretical models are 
also compared to this new treatment of the writhe, and a value of the torsional 
modulus is computed. 

I will use the notations: f3 = l/fceT, F = fk-^T the force exerted on the 
ends of the molecule and 9 the angle of a vector with the direction of the force, 
L the length of the molecule, A = ipk^T the bending modulus, C — itkBT the 
torsional modulus, T = ^k^T the torque exerted on the free end of the molecule, 
z its vertical extension and x the rotation angle of the free end. The torsion 
u){s) of the DNA molecule is the difference per unit length of the rotation angle 
of one strand around the other with respect to the unconstrained state in the 
same conformation. It is a function of the arc length s and can be seen as "twist 
angle density" : The torsion integrated along the molecule is called twist angle 
and is noted fi. 

2 DNA elasticity models 

Thanks to its double helix structure, DNA is a very stable polymer. This sta- 
bility is needed to conserve genetic materials. It gives opportunity to submit 
a single molecule to mechanical constraints (forces above 0.04 pN) without de- 
stroying it and to measure its response. The extension z of a single molecule 
submitted to a force F = fk-eT is the first quantity that has been studied. It 
has to be interpreted within a polymer elasticity model. 

The freely jointed chain (FJC) describes a polymer as a succession of inde- 
pendent sticks of length h. For one stick the elementary partition function is 

3 = y'^cxp(-/6cos0) d(cos0) = (1) 

For a whole chain, the free energy is thus J^fjc — -^sllt in^, leading to a relative 
extension z/L = fb/3 when / tends to zero, and z/L = 1 — {fb)'^^ for large 
forces. 

The worm-like chain is on contrary a continuous model based on the bending 
energy of the axis. The state variable is the tangent unit vector, described by 
its angle 9 with the force, thus the phase space is made of functions 9{s). The 
resistance of the molecule to bending creates a correlation between tangent 
vectors along a curvilinear distance ip ~ 50 nm called the persistence length. 
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The Hamiltonian for this model is obtained by analogy with a quantum system 
(see Ref. |7]): 

/3Hwo™(/^p) = -^|sin^|-/£peos^. (2) 

When / tends to zero the worm- like chain extension is z/L — |/^p. In the large 
force approximation, the molecule is almost aligned with the force (6* <C 1) and 
the Hamiltonian Q reduces to an harmonic oscillator whose ground state gives 
the free energy |S] 

/3^worm = (-flp + y^) ^ [ftp » 1). (3) 



^ V J^p; ~ w^p 

From (3dJ^v/oim/d{f£p) — ~z/ip, when f£p 1, one deduces 

^ (ftp » !)• (4) 



I now compare the two models. At low force, entropy dominates and the 
models both describe an object following Hooke's law. Requiring the relation 
b = 2£p between their parameters makes them equivalent. This is known to be 
a good approximation when f£p = fb/2 < 1 or in other terms F < 0.08 pN. On 
the contrary at large force, the two models do not have the same asymptotic 
behaviour. Comparison between experimental data and theoretical models indi- 
cate that the worm-like chain provides a more appropriate description of DNA 
elasticity 

The double helix structure of DNA has also for consequence a resistance to 
torsion along its axis. Mezard and Bouchiat 7 showed that taking into account 
the local torsion of DNA introduces new elastic terms in the Hamiltonian (O, 
proportional to 7^: A term due to local torsion energy, /3Htorsion = ^^7^ ^74, 
and a term related to the geometry of the molecule 

/?Hrod(/^p, 7) = /3Hworm(/^p) + /^Hfrsion - ^ ^ • (5) 

Z i + COS u 

This improved model is called the rod-like chain. In the large force approxima- 
tion, the denominator 1 -I- cos 6* ~ 2, thus to second order in 6 the Hamiltonian 
rewrites /3Hrod(/^p, 5) ^ /3Hworm(/^p - 7^/4) + /SHtorsion - 7^/4- This approx- 
imation is valid when 9 is small, thus in a domain where the relation Q is 
correct. One deduces the expression of the free energy of the rod-like chain in 
the large force regime: 

/3^rod(/^p, 7) = /3^won„(/^p ~ 7V4) - + (6) 

3 The writhe and its fluctuations 

Experimental devices also allow to change the angle Xj ^-nd to measure how 
the molecule responds to this constraint 0. One says that the molecule is 
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supercoiled. Supercoiling is quantified by the writhing angle of the DNA axis C 

uniini: 

$ = i (/ dr • ® dr' X — — (7) 

2/c Jc \\r-r'f 

This expression is a priori valid only for closed chains. However, in experimental 
devices the molecule can not go around its ends 0] , which is equivalent to having 
a very long molecule only manipulated in a small region. The molecule axis can 
therefore be imaginarily closed and the formula iQ extends to open chains in 
experimental conditions Therefore, the fundamental geometrical relation 

between x, ^ and $ [T^ . 

X = + (8) 
is exact in experimental conditions. Formula © is widely known by DNA 
specialists, in biology and physics ^U]. It explains that a DNA molecule has 
two ways to deal with an applied torque : to modify its local torsion (modify 
J7) or to change its shape (modify 

The twist angle ft is only related to the torsion. The writhe angle $, as 
it is computed from the shape of the axis, is only related to the curvature. A 
general study of DNA elasticity performed by Marko and Siggia asserts that a 
coupling between bending and torsion is at least of the third order in strains, and 
the estimated coupling constant is small . As one wants to focus on regimes 
close to relaxed state, one shall neglect this coupling and assume that there are 
no correlations between fi and (The same assumption has been made in 
Ref. [El CI-) Eq. then gives the fluctuations of x hi this approximation: 

ix') = (n^) + {^') . (9) 

The fluctuations of fi, thanks to the locality of torsional energy and indepen- 
dance of Htorsion with 6, are proportional to the length L, (f^^) = L/£f If 
one measures the torsional modulus of a molecule only taking the angle x into 
account, one one obtains an effective value, denoted k^Tlt^^, related to the real 
torsional constant, fceT^t, through the fluctuations of the writhing angle: 



L L 



(10) 

The effect of writhing of the molecule is then to reduce the measured vale of 
the torsional modulus between its ends. To deduce the value of C along the 
molecule axis from measurements, one needs to know the writhe fluctuations. 

For the rod-like chain model developed in the preceding section, one remarks 
that (x^) is deduced from J^rod, Eq. ©, by differentiation —2f3dJ^i-od/d{'j'^)\j=o = 
(x^), so 

(*'>^i(i-j)^ (11) 

This result, combined with relation was first given by Moroz and Nelson |S] 
as a correction like in Eq. (|10|l . when the force is large enough: 

(^^)=^—^. (12) 
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Their model is therefore an approximation of the more general rod-like chain 
model described by Mezard and Bouchiat. Let us now evaluate the minimum 
force on their common validity range. Expression Hll|) is obtained when the 
molecule is almost straight, namely when f£p ^ 1. That might be translated 
into physical units to F ^ 0.1 pN, thus for forces of the order of 1 pN. 

While Mezard and Bouchiat's approach is more general, it has been pointed 
out that it suffers from a pathology related to the writhe formulation More 
precisely, the expression for the writhe used in Ref. 0, 

J{l-cos9{s))d(f>{s), (13) 

is equal to formula (TJ modulo Att (0 is the azimuthal angle of the tangent vec- 
tor). They are equal only when the molecule can be straightened out without 
cutting it, nor having a point passing through = tt at any time. Otherwise 
the formulae differ by a multiple of 47r ^| . When the molecule is deformed 
into a blob, such configurations where the formulae disagree are numerous. Typ- 
ically, in this case the molecule has loops. The following section is dedicated to 
an estimation of the number of loops in order to extend the validity range of 
formula Hll|l . 

4 Estimation of the loops contribution 

Each loop contributes to the writhing angle <i> by an amount of the order of one 
turn. I will note the mean square value of a loop contribution. This quantity 
will be numerically estimated in the next section. Under the assumption that 
loops are uncorrelated, the central limit theorem asserts that the contribution of 
loops to (<i>^) is n(/£p)A^, where n{fip) is the average number of loops for f£p 
fixed. In Ref. a treatment of loops was proposed, but only equilibrium 
loops were considered. This approach does not explicitly provide any loop size. 
Since the length of one stick being 2£p, the model avoids the smaller loops, that 
require much bending energy. 

The average number of loops is estimated here by defining a loop as a region 
where cos 6* < 0, which gives, using Eq. ^ and b = 2£p, 

nifip) = ^ - r exp(-/6cos0) d(cos0) ~ Aexp(-/£p). (14) 

Thus one obtains an estimate for ($'^) in an extended force range by adding 
the contribution of the loops (regions where cos 9 < 0), given by the preceding 
formula, to the one of the other regions (where it is supposed that 6* ^ 1), given 
by the expression Hll() : 

($2) ^ i (l - 1) ^ + Ia^ exp(-/£p)^ (/^p < 1). (15) 

Adding these contributions is allowed here thanks to Fuller's formula for writhe 
change during a deformation with fixed boundary conditions (see Ref.|14p. 
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When f£p is small, the last term of Eq. (|15|l becomes dominant: The writhe is 
dominated by the loops contribution. is computed in the next section. 

5 Numerical results 

A Monte-Carlo simulation of semiflexible chains was performed (described in 
details in Ref. and the writhe angle of each configuration of the simulated 
equilibrated ensemble was computed, with both expressions and (fU^ . It 
has been shown numerically that extra closure terms are negligible jjl] . Results 
are displayed in figure ^ The estimate of Eq. (fT5|l is also displayed for the 
best value of = 6.4 (one could conjecture that the geometric exact value be 
2tt). The validity ranges estimated above appear reasonable, and the domains 
where different model agree are clearly observed. The distribution of writhe 
angle is Gaussian, of width proportional to L/ip and centered around zero. 
The estimated loops contribution formula agrees quite well with the results 
and explains the difference to other models: It is not counted in one case and 
overestimated in the other one jJH]. If one removes the curves that have at 
least one backfacing region, which corresponds to the definition of a loop (see 
Eq. id}), one gets a curve similar to Moroz and Nelson's. The three models 
agree for F > 1 pN, so the origin of model disagreement is located at low force. 

The numerical results for the writhe angle displayed in figure ^ combined 
with the experimental data of references [3 ^ obtained at low force, in the 
elastic regime, give the value 

h = C/k^T = 93± lOnm. 

This value has to be compared with the ones given in the introduction. The 
theory of Ref. 7^ was applied for forces below this value, in a domain where it 
overestimates the writhe fluctuations. The value obtained by Moroz and Nelson 
is the closest to the recent value of Ref. 0. It was obtained by eliminating data 
points obtained with a force lower than 0.3 pN therefore in a domain where all 
models are equivalent. It is suggested in Ref. that applying a large force 
could lead to an overestimated value for C because of structural modifications 
of the double helix. This question is still open until now. 

6 Discussion 

Unlike Mezard and Bouchiat's theory, neither Moroz and Nelson's model nor 
this work needs to be regularized by any cut-off. The value of this cut-off is 
difficult to relate to an independently measurable quantity. It is suggested in 
Ref. |S] that a treatment of self-avoidance could be necessary. 

In the low force regime, where the freely jointed chain model is valid for DNA, 
let us consider one molecule as a chain of sticks of length h — 2£p and diameter 
d = 2nm. The molecule is located in a region of size R ^ \/Lb. In the low 
force regime, the sticks have almost random orientations, so excluded volume 
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interactions can be estimated with the second virial coefficient. It follows that 
excluded volume interactions appear when p^R^lP'd > 1 where p = {L/b)R~^ is 
the density of sticks. This gives L/b > {b/d)'^ = 2500. The length of the widely 
used A-phage DNA is 16 /xm ~ 300^p|31. For larger forces, the molecule gets 
an orientation and self-avoidance effects are still smaller. It is established in 
Ref. |S| that self-avoidance effects can be neglected as long as z/L> 0.25. This 
argument does not apply as is in the present study since DNA stores torsion. It 
is known, for example in the case of a plectoneme, that self-avoidance stabilizes 
some torque constrained configurations. I focused on the elastic regime, where 
the torque is small and in presence of a small but non-zero force, then applying 
a torque diminishes z/ L. In this situation, the assumption that has been made 
consists in considering self-avoidance negligible when 7 — > 0, or quantitatively 
when z/L > 0.25. As a consequence of self-avoidance the statistical weight of 
conformations with large writhe increases, in other words is underestimated 
in our model. Consequently, taking into account self-avoidance effects would 
result in an small increase of the value of C given above. 

In the numerical work, self-avoidance was not taken into account, following 
the considerations of the preceding paragraph. As a consequence, the studied 
ensemble contains knotted configurations. I showed from the same numerical 
simulations that a 8^p long worm-like chain has a probability of (5 ± 1) 10~^ 
to knot ^Jj. For a force of 0.02 pN the polymer statistic can be split into 
independent elastic blobs of size 8^p. Then the probability to have a knot in at 
least one blob is around 2.10~^ for a 16 fim DNA molecule. If the force is higher, 
the knot probability decreases, as blobs become smaller, and is negligible in the 
experimental force range. Knotted configurations are then believed to play no 
role in the experimental force range. 

In this work, I have investigated the connections between different DNA 
models and showed that Moroz and Nelson's model and Mezard and Bouchiat's 
are related in a simple way, which depends on how DNA loops are taken into 
account. The estimates of the writhe fluctuations of the molecule in those models 
have been compared to numerical results. I have shown that the differences 
between those models are due to the absence (in Moroz and Nelson's model) 
or the overestimate (in Mezard and Bouchiat's model) of the contribution of 
the molecule loops to the writhe. This model takes more accurately writhe 
fluctuations into account whithout introducing any cut-off. It also provides a 
reasonable value for the torsional modulus of DNA under small constraints. 

The author would like to thank A. C. Maggs for introduction to this field 
and discussions, and ENS-Lyon for financial help. 
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Figure 1: Writhe angle fluctuations as a function of the tension. The experimen- 
tal force range is indicated by the grey bar. The plain line presents the results 
of the present letter model, compared to modified Moroz and Nelson's result 
(<I>^) i-p/L = (1 — z/L)/2 (dashed Hne) and Mezard and Bouchiat's (dashed- 
dotted line). Moroz and Nelson's result is valid only for large forces [F > 1 pN), 
validity range of Mezard and Bouchiat's is wider but fails for the lowest exper- 
imental force values, around F = 0.2 pN. The estimated contribution of loops 
is displayed as a dotted line (see text). It corresponds to the best fit for 
in Eq. H15(l |16| . The grey line is obtained from the same data set, but after 
removal of all configurations that have backward facing regions ; therefore there 
are no loops in this restricted ensemble. The results of Mezard and Bouchiat's 
model and the one presented in this letter have been obtained numerically, from 
equilibrated sets of 50.000 chains of 32 ip for each point, with discretisation 
^p/30, error bars are thiner than the line thickness. The two other curves are 
analytical expressions. 
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